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EARLY-PHASE SPECTRA OF "HYPERNOVA" SN 2002AP 
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ABSTRACT 

The spectral evolution of the peculiar SN Ic 2002ap during the first 40 days is presented. The spectra 
display very broad absorption features, which are typical of " hypernovae" . The maximum expansion 
velocity measured on the earliest spectra exceeds 3 x 10 4 km s _1 . The spectrum of SN 2002ap at the 
epoch of maximum brightness resembles that of SN 1997ef more than that of SN 1998bw. The spectral 
evolution of SN 2002ap proceeds at about 1.5 times the rate of SN 1997ef. The parameterized supernova 
spectrum synthesis code SYNOW was used to perform line identification and deduce velocity information 
from the early-phase spectra, which are heavily affected by line blending. The photosphcric velocity, 
as deduced from the fitting results and from the blucshift of the Si II A6355 absorption minimum, is 
lower than in previously studied hypernovae. At advanced epochs, the Si II A6355 absorption minimum 
becomes difficult to distinguish. This may be caused by the growth of [O I] AA6300, 6364 emission. 
Together with the rapid spectral evolution, this suggests that SN 2002ap should enter the nebular phase 
sooner than previously studied hypernovae. 



Subject headings: gamma rays:bursts — line:identification- 
— supernovae:individual(SN 2002ap) 



-supernovae : general 



1. INTRODUCTION 

The very powerful supernova (SN) explosions called "hy- 
pernovae" (Iwamoto et al. 1998) are perhaps the most 
energetic events in the current universe. Their kinetic en- 
ergy is ~ 10 52 ergs, which is about 10 times larger than of 
normal SNe. SN 1998bw, probably associated with GRB 
980425 (Galama et al. 1998), was the first recognized hy- 
pernova. Although classified as a Type Ic SN (SN Ic), the 
spectrum of SN 1998bw showed unusual broad spectral 
features in the early phase. The derived velocities were as 
large as 30,000 km s" 1 . The kinetic energy of SN 1998bw 
was deduced to be 3 - 5 x 10 52 ergs (Iwamoto et al. 1998; 
Nakamura et al. 2001). Other Type Ic hypernovae were 
subsequently recognized based on their broad absorption 
features. These include SNe 1997ef (e.g., Iwamoto et al. 
2000), 1998ey (Garnavich, Jha, & Kirshner 1998), and 
the recently discovered 2002bl (Filippenko et al. 2002). 
Moreover, SN 1997dq was suggested to be a hypernova 
from its similarity to SN 1997ef at later phases (Matheson 
et al. 2001). The kinetic energy deduced for SN 1997ef 
was about 1.8 x 10 52 ergs (Mazzali, Iwamoto, & Nomoto 
2000), which was smaller than that of SN 1998bw, but still 
significantly larger than normal SNe. Hypernovae also oc- 
cur among SNe Iln, for example, SN 1997cy, one of the 
brightest SNe ever discovered and characterized by a large 
kinetic energy, ~ 3 x 10 52 ergs (Turatto et al. 2000), and 
SN 1999E (Cappellaro, Turatto, & Mazzali 1999). The 
range of properties of hypernovae, and their possible asso- 
ciation with GRBs are also topics of great interest (e.g., 
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Nomoto et al. 2001). 

SN 2002ap was discovered in M74 on January 29.4 UT 
by Y.Hirose (Nakano et al. 2002). Spectra obtained on 
January 30 - 31 showed very broad absorption features, 
resembling those in SNe 1998bw and 1997ef. Therefore, 
SN 2002ap was suggested to be a hypernova (Kinugasa et 
al. 2002; Meikle et al. 2002; Gal- Yam & Shemmer 2002). 
SN 2002ap is the nearest hypernova discovered to date, 
hence it was a good target even for small telescopes. Vari- 
ous multi-wavelength observations of SN 2002ap have been 
performed (summarized in Gal- Yam, Ofek, & Shemmer 
2002). NIR photometry was also performed frequently at 
Gunma Astronomical Observatory (GAO) (E. Nishihara 
et al. 2002, in preparation). 

We report here on the optical spectroscopy of SN 2002ap 
performed at GAO. We compare the spectra of SN 2002ap 
with those of other hypernovae. We also used the param- 
eterized SN synthetic-spectrum code SYNOW to perform 
line identification and deduce velocity information from 
the spectra in the early phase, when very broad features 
resulting from line- blending were present. 

2. OBSERVATIONS AND DATA REDUCTION 

Observations of SN 2002ap started on the night of the 
discovery report. Nightly early phase spectral observations 
were carried out at GAO and Bisei Astronomical Obser- 
vatory (BAO) for about 40 days, from discovery to solar 
conjunction. Spectra were obtained with Gunma Compact 
Spectrograph (GCS) attached to the classical-Cassegrain 
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focus of the GAO 65cm telescope (F/12) and with a spec- 
trograph attached to the BAO 1.01-m telescope (F/12). 
The spectral resolving power (A/AA) of GCS is about 500 
at 5500 A and the wavelength coverage is from 3800 to 
7600 A. The BAO spectrograph covers from 4700 to 7000 
A with a resolving power of about 1000. Observational pa- 
rameters are listed in Table 1. Since all the spectra were 
obtained with narrow slits (2" on GCS and 2.6" at BAO, 
compared to a typical seeing of 2-3") and the sky condi- 
tions were generally non-photometric, the flux was not on 
an absolute flux scale. The slit was aligned in E-W direc- 
tion at GAO and in the N-S direction at BAO, which in 
both cases was not the parallactic angle. Consequently, 
the blue end of the spectra could be significantly affected 
at advance epochs, as the observations were done at low 
altitude. In those cases, we cut the blue end of the spectra 
(see Figure 1). 

Data reduction was performed with IRAF 7 . Calibration 
for the relative intensity was performed based on the obser- 
vations of spectrophotometric standard stars. A Fe-Ne-Ar 
lamp for GAO and a Fe-Ne lamp for BAO were used for 
the wavelength calibration. 

EDITOR: PLACE TABLE 1 HERE. 



3. RESULTS AND DISCUSSION 

The time-sequence of the calibrated spectra of 
SN 2002ap is shown in Figure 1. The evolution of the 
spectra is apparent. As time progresses, most absorption 
features shift to the red, and the continuum becomes red- 
der. 

EDITOR: PLACE FIGURE 1 HERE. 

The first two spectra, taken on the same night, show a 
blue continuum, with neither deep absorptions nor strong 
emissions (see also the top of Figure 2). Broad and shal- 
low depressions are seen near 4700, 5700, and 6200 A. 
These broad features resemble those of SNe 1998bw and 
1997ef in the earliest phases (Patat et al. 2001; Iwamoto 
et al. 2000). However, SN 2002ap is bluer, possibly be- 
cause it was caught at an earlier epoch. The fact that the 
first observations must have been very close to the time 
of explosion is exemplified by the rapid development of 
spectral features. The first two spectra also show a broad 
peak near 5000A. But in the third spectrum taken the fol- 
lowing night, this peak has shifted significantly to the red 
(~ 5100A), indicating a rapid decrease of the material ve- 
locity sampled by the spectra, as the photosphere moves 
inwards. 

Photometric observation at Wise Observatory estab- 
lished that the peak of the light curve occurred on Febru- 
ary 7.1 and 8.8 UT in the B and U-band, respectively 
(Gal- Yam et al. 2002). We designate the epochs using 
the number of days from B maximum. Figure 2 also shows 
the spectrum of SN 2002ap taken on day +2 (February 9.4 
UT), near the time of maximum. We compared this with 
the near-maximum spectra of other SNe Ic: SNe 1998bw, 

7 IRAF is distributed by the National Optical Astronomy Observato 
Astronomy, Inc., under cooperative agreement with the National Sci 



1997ef (hypernovae) and 19941 (a normal SN Ic; Nomoto et 
al. 1994). Based on the line width, the spectra can easily 
be separated into two groups. One is that of "hypernovae" , 
including SNe 1998bw, 1997ef, and 2002ap, while the other 
includes only the normal Type Ic SN 19941. Therefore, 
SN 2002ap is easily identified as a hypernova. Moreover, 
the spectrum of SN 2002ap shows emission-like features 
around 4600 A and 5300 A, resembling the spectrum of 
SN 1997ef more than that of SN 1998bw. 

EDITOR: PLACE FIGURE 2 HERE. 

In Figure 3 we compare similar-looking spectra of 
SNe 2002ap and 1997ef. We can see that the spectral 
evolution of SN 2002ap is about 1.5 times faster. The 
earliest available spectrum of SN 1997ef, on day —6, ap- 
pears to have properties intermediate between the spec- 
trum of SN 2002ap on day —1 and that on day —6, judg- 
ing from the progressive redshifting of line features. The 
near-maximum spectrum of SN 1997cf resembles that of 
SN 2002ap on day +2 in all major features. The simi- 
larity between the spectrum of SN 1997ef on day +19 and 
that of SN 2002ap on day +15 is striking, both showing mi- 
nor features of comparable intensities at about 4900, 5900, 
6200, 7100 and 7300 A. Most of the new features visible in 
the spectrum of SN 1997ef on day +52, some of which are 
probably net emissions, are already visible in SN 2002ap 
on day +25, although they are not as developed. Both of 
these spectra suggest that the SNe are making a transition 
to the nebular phase, while the spectrum of SN 1997ef on 
day +27 is still essentially photospheric in nature (Mazzali 
et al. 2000). 

EDITOR: PLACE FIGURE 3 HERE. 

Especially in the earliest spectra, most features are so 
broad because of line blending that it is not easy to iden- 
tify which lines contribute to them. Thus we used the 
parameterized SN spectrum synthesis code SYNOW, de- 
scribed by Fisher (2000), to perform line identification 
and deduce velocity information. It is simpler and faster 
than the Monte-Carlo code of Mazzali et al. (2002), and 
better suited for our objective, as it does not require the 
construction of realistic SN models. In this paper, the 
radial-dependent line optical depth above the photosphere 
is assumed to be proportional to (v/f p h)~", where v is 
the expansion velocity and u p h is the photospheric veloc- 
ity. A reasonable fit to the observations could be obtained 
using just a few ions: Ca II, O I, Si II, Fe II, Co II and 
Ni II, consistent with the line identification in SNe 1997ef 
and 1998bw (Branch 2001). Synthetic spectra and line 
identification for day —7 and +2 are shown in Figure 2. 

Either a high w p h or a small n can make the lines broad. 
Guided by Branch's (2001) work, we found that both are 
necessary to account for the extreme line widths. Similar 
results were obtained by Mazzali et al. (2002), who used 
a Monte-Carlo code to synthesize some selected spectra of 
SN 2002ap. Our best fit value of the power-law index was 
3 for the first spectrum. This index value was used for 
the radial dependence of the opacity for all lines and at all 
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epochs. As shown in Figure 4, the photosphcric velocity 
obtained from the fitting decrease from 35,000 km s" 1 on 
day —7 to 13,000 km s _1 on day +2, with the uncertainty 
being about 1000 km s _1 empirically (Branch et al. 2002). 

Our fits are reasonable for these epochs, allowing us to 
identify most absorption lines. However, some problems 
remain. One is the unwanted "peak" around 4100 A, which 
is actually the remainder of continuum. We can remove it 
by introducing both Ti II and Sc II lines. However, there 
is no other argument for Sc II in SNe Ic. Such difficulties 
may show the limitations of SYNOW. They are probably 
due to the assumption of a pure resonant scattering source 
function and on the simple power-law dependence of the 
optical depth on velocity applied to all ions. The hump 
around 5200 A is also incorrectly reproduced, probably 
the consequence of not transferring sufficient flux from the 
blue (~ 4000A) in SYNOW. Finally, we introduced Na I 
to fit the absorption around 5700 A. He I A5876 may also 
contributes there. The NIR spectra show a significant fea- 
ture, possibly due to He I A10830 (Motohara et al. 2002). 
However, that feature in SN 19941 might have been C I, 
O I, and/or Si I lines. (Baron et al. 1999; Millard et al. 
1999). 

Among those absorption lines identified, the Si II 
AA6347, 6371 (collectively called A6355) absorption min- 
imum was easily distinguishable during almost the entire 
observation period. Therefore, we measured the minimum 
of the Si II absorption as an alternative method of estimat- 
ing the photosphcric velocity. In Figure 1, we marked with 
an arrow the minimum of the Si II absorption. The evo- 
lution of the photospheric velocity, as deduced from this 
method, is also shown in Figure 4. At the earliest phases, 
the typical statistical error was around ± 3000 km s _1 . 
Moreover, strong line blending makes the estimate of the 
velocity somewhat uncertain at those epochs. Near max- 
imum light, although the minima arc clearly visible (the 
statistical error was about ± 1500 km s _1 ), the system- 
atic error on the velocity is mostly due to the ambiguity 
introduced by the uncertain continuum. At late phases, 
the error (around ± 2500 km s _1 ) was mainly due to poor 
statistics. The difference between the velocities from Si II 
and those from SYNOW fitting may be mostly caused by 
systematic errors introduced by the uncertainty in the con- 
tinuum and the heavy line blending. 

EDITOR: PLACE FIGURE 4 HERE. 



In Figure 4, we compare the evolution of the photo- 
spheric velocity of SN 2002ap, 1998bw and 1997ef (Patat 
et al. 2001). We adopted for SN 2002ap an explosion date 
of January 28.9 UT (day —9.2), as proposed by Mazzali 
etal. (2002). The photospheric velocity of SN 2002ap de- 
clines faster than for the other two hypernovae, suggesting 
that the ejecta mass of SN 2002ap is smaller. Although the 
velocity evolution of the three SNe is initially well fitted 
with exponential functions, both SNe 1997ef and 1998bw 
eventually seem to flatten out, SN 1998bw at about 7000 
km s- 1 and SN 1997ef at about 2000 km s" 1 . While the 
measurement for SN 1998bw can be done directly on the 
spectra, the value for SN 1997ef is derived from model 
fitting, since, as for SN 2002ap, the Si II line becomes 
indistinguishable at advanced epochs (see spectra on day 
+52 for SN 1997ef and on day +25 for SN 2002ap in Fig- 
ure 3). No sign of the flattening for SN 2002ap was seen 
over the period of our observations, but it might occur 
later, which may have to be established through detailed 
spectral modeling. 

The rapid evolution of the photosphcric velocity and of 
the spectra of SN 2002ap suggests that SN 2002ap will 
develop a nebular spectrum earlier than the other hyper- 
novae presented here. In SN 1997ef nebular emission was 
observed in the Ca II IR triplet on day +27. Our spec- 
tra of SN 2002ap unfortunately do not extend to that re- 
gion. However, the spectrum on day +25, which is inter- 
mediate between the day +27 and the day +52 spectra 
of SN 1997ef, may show' the onset of weak [O I] AA6300, 
6364 net emission. E. Nishihara et al. (2002, in prepara- 
tion) show that a break occurs in the evolution of the IR 
colors around February 25, suggesting the onset of nebu- 
lar emission. The first spectra of SN 2002ap after solar 
conjunction, taken on day 131 and 140, show that it has 
already entered the nebular phase (Leonard et al. 2002). 

The broad spectral features, suggesting high kinetic en- 
ergy, distinguish hypernovae from normal SNe Ic. How- 
ever, the observations of SN 2002ap suggest that there is a 
variety of hypernovae, based on the observational features 
(spectra, spectral evolution, luminosity, light curves, the 
association with a GRB, etc.). More observations of hyper- 
novae are necessary in order to understand what physical 
parameters are responsible for the variations among hy- 
pernovae and whether there is a gap in properties between 
hypernovae and normal SNe Ic (Figure 2 of Nomoto et al. 
2002). 
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Fig. 1. — Spectral evolution of SN 2002ap from -7 to +30 days relative to B-maximum (Feb. 7.1; Gal- Yam ct al. 2002). The spectra have 
been shifted vertically by arbitrary amounts. The horizontal axis is the observed wavelength. Arrows mark the absorption minima of Si II 
A6355 (See text on the errors). The dashed line marks the rest wavelength (at the redshift of M74) of Si II. The blue ends of advanced-epoch 
spectra, which are affected by observing at low altitude, have been cut for clarity. The symbol © marks the telluric B-band absorption. 
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Fig. 2. — The spectra of SN 2002ap taken on Jan. 31 and Feb. 9 are compared with a synthetic spectrum computed with SYNOW, and the 
line identification is shown. Also shown are the maximum-light spectra of other SNe Ic for comparison (SNc 1998bw and 1997cf from Fig.l 
of Mazzali et al. 2002 and SN 19941 from Fig. lb of Millard et al. 1999). The spectra have been scaled and shifted vertically by arbitrary 
amounts, and shifted in wavelength to the rest-frame of the parent galaxies. The observed spectra of SN 2002ap have been corrected for the 
interstellar extinction, E(B — V) = 0.09 (Takada-Hidai, Aoki, & Zhao 2002) and for the host galaxy recession, 657 km s _1 . 
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Fig. 3. — A comparison of various spectra of SN 2002ap and similar-looking spectra of SN 1997ef. The spectrum of SN 1997ef on day ~ 
—6 is inserted between two adjacent SN 2002ap spectra. Other spectra of SN 1997ef are polluted by sharp H II emissions and hence easily 
recognizable. The spectra have been scaled and shifted vertically by arbitrary amounts, and shifted in wavelength to the rest-frame of the 
parent galaxies. 



8 



Kinugasa et al. 



O 



• i — i 
o 
o 

> 



o 



o 
o 
o 




O SN 2QQ2ap (from Sill) 

• SN 2002ap (from fitting) 

□ SN 1997ef (from Sill) 

■ SN 1997ef (from fitting) 

SN 1998bw (from Sill) 







To 

□ 



<> <> 







o 



□ 



- - - - < 


> < 


> 









o 

00 



o 



20 



40 



60 



Days From Explosion 



Fig. 4. — The evolution of the photosphcric velocity of hypcrnovae. The SN 2002ap data arc deduced from the Si II A 6355 absorption 
minima (Fig.l) and from SYNOW fitting. The explosion date was set at Jan. 28.9 (Mazzali ct al. 2002). For comparison, the data for the 
other hypernovae are shown from Patat et al. (2001) for SN 1998bw, and Mazzali ct al. (2000) for SN 1997ef. The curves arc fitted with 
exponential functions for SNe 2002ap and 1997ef (Si II data), and with exponential functions plus constants for SNe 1998bw and 1997cf 
(fitting data). 
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Table 1 
Observation log. 



Date (UT) a 


Phase (Days) b 


Exposure (s) 


Tclcscope c 


2002 Jan. 


31.40 


-7 


3600 


G 


2002 Jan. 


31.47 


-7 


900 


B 


2002 Fob. 


1.40 


-6 


3600 


G 


2002 Feb. 


6.39 


-1 


5400 


G 


2002 Feb. 


6.41 


-1 


2700 


B 


2002 Feb. 


7.41 





3600 


G 


2002 Feb. 


9.39 


+2 


3600 


G 


2002 Feb. 


10.41 


+3 


2400 


G 


2002 Feb. 


12.39 


+5 


3600 


G 


2002 Feb. 


13.40 


+6 


2400 


G 


2002 Feb. 


14.40 


+7 


3600 


G 


2002 Feb. 


16.39 


+9 


3600 


G 


2002 Feb. 


20.39 


+ 13 


3600 


G 


2002 Feb. 


21.42 


+ 14 


2400 


G 


2002 Feb. 


22.40 


+ 15 


3600 


G 


2002 Feb. 


25.42 


+ 18 


2400 


G 


2002 Mar. 


4.41 


+25 


2400 


G 


2002 Mar. 


8.40 


+29 


2400 


G 


2002 Mar. 


9.41 


+30 


1800 


G 





a Starting time of the exposures. 

b Days from B-maximum : Feb. 7.1 UT (MJD=52312.1) (Gal- Yam 
ct al. 2002). This was about 9.2 days after the explosion, which 
estimated by Mazzali et al. (2002) to have occurred on Jan. 28.9. 

C 'G' and 'B' represent GAO 65cm and BAO 1.01m telescopes, 
respectively. 



